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Abstract. The main objective of this article is to derive a new set of gravi- 
tational field equations and to establish a new unified theory for dark energy 
and dark matter. The new gravitational field equations with scalar potential 
tp are derived using the Einstein-Hilbert functional, and the scalar potential 
tp is a natural outcome of the divergence- free constraint of the variational ele- 
ments. Gravitation is now described by the Riemannian metric gij, the scalar 
potential tp and their interactions, unified by the new gravitational field equa- 
tions. Associated with the scalar potential tp is the scalar potential energy 

4 4.. 

density $ = g C w Q 9 tJ DiDjtp, which represents a new type of energy caused 
by the non-uniform distribution of matter in the universe. The negative part 
of this potential energy density produces attraction, and the positive part pro- 
duces repelling force. This potential energy density is conserved with mean 

4 

zero: J M &dM = 0. The sum of this new potential energy density ^ G <E> 
and the coupling energy between the energy-momentum tensor T;j and the 
scalar potential field tp gives rise to a new unified theory for dark matter and 
dark energy: The negative part of this sum represents the dark matter, which 
produces attraction, and the positive part represents the dark energy, which 
drives the acceleration of expanding galaxies. In addition, the scalar curvature 
of space-time obeys R = 8 ^ 4 G T + <!>. Furthermore, the new field equations 
resolve a few difficulties encountered by the classical Einstein field equations. 
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1. Introduction and Summary 

The main aim of this article is an attempt to derive a new theory for dark matter 
and dark energy, and to derive a new set of gravitational field equations. 

The primary motivation of this study is the great mystery of the dark matter and 
dark energy. The natural starting point for this study is to fundamentally examine 
the Einstein field equations, given as follows: 

1 8ttG 

(1.1) R ij --g ij R~ -Tij, 

where Rij is the Ricci curvature tensor, R is the scalar curvature, gij is the Riemann- 
ian metric of the space-time, and Ty is the energy-momentum tensor of matter; 
see among many others [T]. The Einstein equations can also be derived using the 
Principle of Lagrangian Dynamics to the Einstcin-Hilbert functional: 

(1.2) F( 9ij ) = [ R^dx, 

J M 

whose Euler-Lagrangian is exactly Rij — ^gijR, which is the left hand side of the 
Einstein field equations (|1.1|) . It is postulated that this Euler-Lagrangian is bal- 
anced by the symmetric energy-momentum tensor of matter, Tj,-, leading to the 
Einstein field equations (jl.ip . The Bianchi identity implies that the left hand side 
of the Einstein equations is divergence-free, and it is then postulated and widely 
accepted that the energy-momentum tensor of matter is divergence-free as well. 
However, there are a number of difficulties for the Einstein field equations: 
First, the Einstein field equations failed to explain the dark matter and dark 
energy, and the equations are inconsistent with the accelerating expansion of the 
galaxies. In spite of many attempts to modify the Einstein gravitational field equa- 
tion to derive a consistent theory for the dark energy, the mystery remains. 
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Second, we can prove that there is no solution for the Einstein field equations 
for the spherically symmetric case with cosmic microwave background (CMB). One 
needs clearly to resolve this inconsistency caused by the non-existence of solutions. 

Third, from the Einstein equations it is clear that 

(1.3) *=^T, 

where T — g lJ Tij is the energy- momentum density. A direct consequence of this 
formula is that the discontinuities of T give rise to the same discontinuities of the 
curvature and the discontinuities of space-time. This is certainly an inconsistency 
which needs to be resolved. 

Fourth, it has been observed that the universe is highly non-homogeneous as 
indicated by e.g. the "Great Walls", filaments and voids. However, the Einstein 
equations do not appear to offer a good explanation of this inhomogeneity. 

These observations strongly suggest that further fundamental level examinations 
of the Einstein equations are inevitably necessary. It is clear that any modification 
of the Einstein field equations should obey three basic principles: 

• the principle of equivalence, 

• the principle of general relativity, and 

• the principle of Lagrangian dynamics. 

The first two principles tell us that the spatial and temporal world is a 4-dimensional 
Riemannian manifold (M, g^ ) , where the metric {g^ } represents gravitational po- 
tential, and the third principle determines that the Riemannian metric {gij} is an 
extremum point of the Lagrangian action. There is no doubt that the most natu- 
ral Lagrangian in this case is the Einstein-Hilbcrt functional as explained in many 
classical texts of general relativity. 

The key observation for our study is a well-known fact that the Riemannian 
metric is divergence-free. This suggests two important postulates for deriving a 
new set of gravitational field equations: 

• The energy-momentum tensor Ty of matter need not to be divergence-free 
due to the presence of dark energy and dark matter; and 

• The field equations obey the Euler-Lagrange equation of the Einstein- 
Hilbert functional under the natural divergence-free constraint, with di- 
vergence defined at the extremum Riemannian metric g: 

DgXki = 0, Xki — Xik- 

Here D l g is the contra- variant derivative with respect to the extremum point 
<?, and Xij are the variational elements. Namely, for any X — {X^} with 



D^Xn = 0, 



f 

x^b A 



lim -\F{ 9ij + XXij ) - F( 9ij )] = (SF( gij ),X) = 0. 



For this purpose, an important part of this article is to drive an orthogonal 
decomposition theorem of tensors on Riemannian manifolds, which we shall explain 
further in the last part of this Introduction. 

Under these two postulates, using the orthogonal decomposition theorem of ten- 
sors, we derive the following new set of gravitational field equations with scalar 
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potential: 

1 &7rG 

(1.4) R VJ - -gijR = — — T ij ~ D.Djtp, 

where the scalar function <p : M — > M is called the scalar potential. 

The corresponding conservations of mass, energy and momentum are then re- 
placed by 

87tG 

(1.5) div (DiDjtp + —T l3 ) = 0, 

and the energy-momentum density T = g^Tij and the scalar potential energy 
density g^L$ = D.Djip satisfy 

(1.6) R=^ T + <Z>, 



(1.7) 



M 



The scalar potential energy density gf^j^ has a number of important physical 
properties: 

1 . Gravitation is now described by the Riemannian metric gij , the scalar po- 
tential ip and their interactions, unified by the new gravitational field equa- 
tions (fL~4|) . 

4 

2. This scalar potential energy density gf^^ represents a new type of en- 
ergy/force caused by the non-uniform distribution of matter in the uni- 
verse. This scalar potential energy density varies as the galaxies move and 
matter of the universe redistributes. Like gravity, it affects every part of 
the universe as a field. 

4 

3. This scalar potential energy density & consists of both positive and neg- 
ative energies. The negative part of this potential energy density produces 
attraction, and the positive part produces repelling force. The conserva- 
tion law (|1.7[) amounts to saying that the the total scalar potential energy 
density is conserved. 

4 

4. The sum of this new potential energy density and the coupling energy 
between the energy-momentum tensor T^- and the scalar potential field <p, 
as described e.g. by the second term in the right-hand side of (jl.9p . gives 
rise to a new unified theory for dark matter and dark energy: The negative 
part of e represents the dark matter, which produces attraction, and the 
positive part represents the dark energy, which drives the acceleration of 
expanding galaxies. 

5. The scalar curvature of space-time obeys (|1.6p . Consequently, when there 
is no normal matter present (with T = 0), the curvature R of space-time is 
balanced by R = Therefore, there is no real vacuum in the universe. 

6. The universe with uniform distributed matter leads to identically zero scalar 
potential energy, and is unstable. It is this instability that leads to the 
existence of the dark matter and dark energy, and consequently the high 
non-homogeneity of the universe. 

Hereafter, we further explore a few direct consequences of the above new gravi- 
tational field equations. 
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First, the new field equations are consistent with the spherically symmetric case 
with cosmic microwave background (CMB). Namely, the existence of solutions in 
this case can be proved. 

Second, our new theory suggests that the curvature R is always balanced by $ 
in the entire space-time by (jl.6p . and the space-time is no longer flat. Namely the 
entire space-time is also curved and is filled with dark energy and dark matter. 
In particular, the discontinuities of R induced by the discontinuities of the energy- 
momentum density T, dictated by the Einstein field equations, are no longer present 
thanks to the balance of $. 

Third, this scalar potential energy density should be viewed as the main cause 
for the non-homogeneous distribution of the matter/galaxies in the universe, as the 
dark matter (negative scalar potential energy) attracts and dark energy (positive 
scalar potential energy) repels different galaxies; see (|1.9|) below. 

Fourth, to further explain the dark matter and dark energy phenomena, we 
consider a central matter field with total mass M and radius tq and spherical 
symmetry. With spherical coordinates, the corresponding Riemannian metric must 
be of the following form: 

(1.8) ds 2 = -e u c 2 dt 2 + e v dr 2 + r 2 {d9 2 + sin 2 Qdip 2 ), 

where u — u(r) and v — v(r) are functions of the radial distance. With the new 
field equations, the force exerted on an object with mass m is given by 



(1.9) F = mMG 







Rr~ 









R = G> for r > r . 



where 6 — 2GM/ c 2 , R is the scalar curvature, and ip is the scalar potential. The first 
term is the classical Newton gravitation, the second term is the coupling interaction 
between matter and the scalar potential ip, and the third term is the interaction 

4 

generated by the scalar potential energy density as indicated in (| !.€>[) (R = $ 

for r > ro). In this formula, the negative and positive values of each term represent 
respectively the attracting and repelling forces. It is then clear that the combined 
effect of the second and third terms in the above formula represent the dark matter, 
dark energy and their interactions with normal matter. 

Also, importantly, this formula is a direct representation of the Einstein's equiv- 
alence principle. Namely, the curvature of space-time induces interaction forces 
between matter. 

In addition, one can derive a more detailed version of the above formula: 



(1.10) F = mMG 



-±+(2+ 6 -)er 2 + ?f+ 1 -(2+ S ~)r> I r rRdr 
r \ r I o o \ r 



where e > 0. The conservation law (|1.7p of $ suggests that R behaviors as r~ 2 for 
r sufficiently large. Consequently the second term in the right hand side of (|1.10[) 
must dominate and be positive, indicating the existence of dark energy. 

In fact, the above formula can be further simplified to derive the following ap- 
proximate formula: 



(1.11) F = mMG 

(1.12) k = 4 x 10~ 1S km~ l , h = 10~ 57 km- 



\ - — + k i r 
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Again, in the first term represents the Newton gravitation, the attracting 

second term stands for dark matter and the repelling third term is the dark energy. 

The mathematical part of this article is devoted to a rigorous derivation of the 
new gravitational field equations. 

First, as mentioned earlier, the field equations obey the Euler-Lagrange equation 
of the Einstein-Hilbert functional under the natural divergence- free constraint, with 
divergence defined at the extremum Riemannian metric g: 

(f .13) (SF( gij ), X)=0 V X = {Xij} with D^Xij = 0, X l3 = X }1 . 

As the variational elements X are divergence- free, (j!.13|) does not imply SF(gij) = 
0, which is the classical Einstein equations. In fact, (jl.131) amounts to saying that 
SF(gij) is orthogonal to all divergence-free tensor fields X. 

Hence we need to decompose general tensor fields on Riemannian manifolds into 
divergence-free and gradient parts. For this purpose, an orthogonal decomposition 
theorem is derived in Theorem 13.11 In particular, given an (r, s) tensor field u € 
L 2 (TJM), we have 

(1.14) u = Vtp + v, divw = 0, tj) G H^iXs-i)- 

The gradient part is acting on an (r, s — 1) tensor field ip- 

Second, restricting to a (0, 2) symmetric tensor field u, the gradient part in the 
above decomposition is given by 

VV> with^ = {^i} eH l (T$M). 

Then using symmetry, we show in Theorem 13.21 that this (0, 1) tensor ip can be 
uniquely determined, up to addition to constants, by the gradient of a scalar field 
tfr. 

and consequently we obtain the following decomposition for general symmetric (0, 2) 
tensor fields: 

(1.15) Uij = Vij + DiDjcp, D i v ij = 0, ipeH 2 (M). 

Finally, for the symmetric and divergence free (0, 2) field SF(gij), which is the 
Euler-Lagrangian of the Einstein-Hilbert functional in (11.131) and is orthogonal to 
all divergence-free fields, there is a scalar field ip £ H 2 (M) such that 

SF(gij) = DiDjip, 

which, by adding the energy- momentum tensor T^-, leads to the new gravitational 
field equations (|1.4[) . 

We remark here that the orthogonal decompositions (|1.14[) and (|1.15[) are rem- 
iniscent of the orthogonal decomposition of vectors fields into gradient and diver- 
gence parts, which are crucial for studying incompressible fluid flows; see among 
many others [8] . 

This article is divided into two parts. The physically inclined readers can go 
directly to the physics part after reading this Introduction. 
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Part 1. Mathematics 



2. Preliminaries 

2.1. Sobolev spaces of tensor fields. Let (M, gy) be an n-dimensional Riemann- 
ian manifold with metric (g'-y), and E = T T S M be an (r, s)-tensor bundle on M. A 
mapping it : M —¥ E is called a section of the tensor-bundle £ or a tensor field. In 
a local coordinate system x, a tensor field u can be expressed component-wise as 
follows: 

u = \ u il'.'.T s ( x ) I 1 <h,--- ,is,ji,--- ,Jr <n\ , 

where w^.'.'.'.f' (x) are functions of x € [/. The section u is called C r -tensor field or 
C r -section if its components are C r -functions. 

For any real number 1 < p < oo, let L P (E) be the space of all L p -integrable 
sections of E: 



L p (E) 

equipped with the norm 



u: M -> E 



\u\ p dx < oo 



M 



\U\\lp = 



] p dx 



M 



1 J2K--i:\ Pd x 

M 



VP 



For p = 2, L 2 (E) is a Hilbert space equipped with the inner product 



(2.1) 



fc,. 



det(<7ij), and y/—gdx is the 



where (gij) is Riemannian metric, (g l 
volume element. 

For any positive integer k and any real number 1 < p < oo, we can also define 
the Sobolev spaces W k,p (E) to be the subspace of L P {E) such that all covariant 
derivatives of u up to order k are in L P (E). The norm of W k ' p (E) is always denoted 

by || • \\w*-p- 

As p — 2, the spaces W k,p (E) are Hilbert spaces, and are usually denoted by 
(2.2) H k (E) = W k ' 2 (E) for/c>0, 

equipped with inner product (•, -) H k and norm || • ||#fc. 

2.2. Gradient and divergent operators. Let u : M — > E be an (r, s)-tensor 
field, with the local expression 

(2-3) u ={<;;£}• 

Then the gradient of u is defined as 

(2.4) v, 

where D = (£>i, • • • ,D n ) is the covariant derivative. It is clear that the gradient 
Vu defined by (|2 ,4[) is an (r, s + l)-tensor field: 

Vu : M V S+1 M. 

We define V*m as 



(2.5) 



V*m = {g kl Diu} : M -> T s r+i Af for u as in J23 
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For an (r + 1, s)-tensor field u = {u^...^ ^ r }, the divergence of u is defined by 

(2.6) divu = {A<::.t' jV }- 

Therefore, the divergence div u defined by (|2.6|) is an (r, s)-tensor field. Likewise, 
for an (r, s + l)-tensor field 

the following operator is also called the divergence of u, 

(2.7) divu = {D l u£:f.. is }, 

where D l = g lk D^, which is an (r, s)-tensor field. 

For the gradient operators (|2.4[) - (|2.5[) and the divergent operators (|2.6[) - (I2.7[) . it 
is well known that the following integral formulas hold true; see among others [5] . 

Theorem 2.1. Let (M,gtj) be a closed Riemannian manifold. If u is an (r — l,s)- 
tensor and v is an (r, s) tensor, then we have 

(2.8) (V*u,«) = -(u,divv), 

where V*u is as in 12. 5\) and div v is as in 112. 6\) . the inner product (■, ■) is as 
defined by A2.1]) . If u is an (r, s — \)-tensor and v is an (r, s) tensor, then 

(2.9) (Vu,w) = -(u,divv), 
where Vm is as in \2.1$ and divv is as in \2. 7[ ). 

Remark 2.1. If M is a manifold with boundary dM ^ 0, and u\qm = or 
v\dM — 0, then the formulas (|2.8[) and (|2.9[) still hold true. 

2.3. Acute-angle principle. Let H be a Hilbert space equipped with inner prod- 
uct (■, •) and norm || • ||, and G : H — >• H be a mapping. We say that G is weakly 
continuous if for any sequence {u n } C H weakly converging to uq, i.e. 

u n "o in i?, 

we have 

lim (Gu n ,v) = (Guq,v) Vw S ff. 

n— >oo 

If the operator G is linear and bounded, then G is weakly continuous. The following 
theorem is called the acute-angle principle [5]. 

Theorem 2.2. If a mapping G : H — > H is weakly continuous, and satisfies 

(Gu,u) > a\\u\\ 2 - /3, 
for some constants a, (3 > 0, then for any f g H there is a ua G H such that 

(Gu -f, v )=0 Vw G H. 

3. Orthogonal Decomposition for Tensor Fields 

3.1. Main theorems. The aim of this section is to derive an orthogonal decompo- 
sition for (r, s)-tensor fields with r + s > 1 into divergence-free and gradient parts. 
This decomposition plays a crucial role for the theory of gravitational field, dark 
matter and dark energy developed later in this article. 

Let M be a closed Riemannian manifold, and v G L 2 (TJM) (r + s > 1). We say 
that v is divergence-free, i.e., div v = 0, if for all V-0 G L 2 (T^M), 

(3.1) («,VVO=0. 
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Here ip <E H^TJ^M) or H 1 ('IJ_ 1 M), and (-, •) is the L 2 -inner product defined by 

We remark that if v G /^(TJM) satisfies (|3TT]> . then u is weakly differentiable, 
and divu = in the L 2 -sense. If v G L 2 (T^M) is not differential, then (|3.1[) means 
that u is divergence-free in the distribution sense. 

Theorem 3.1 (Orthogonal Decomposition Theorem). Let M be a closed Riemann- 
ian manifold, and u G L (TJJVf) with r + s > 1. T/ie following assertions hold true: 

(1) TTie tensor field u has the following orthogonal decomposition: 

(3.2) u = Vtp + v, 

where ip £ H l (Tl~ l M) or ip G ff 1 (T a r _ 1 M), and div u = 0. 

(2) // M compact, then u can be orthogonally decomposed into 

(3.3) u = V(p + v + h, 

where ip and v are as in tS.ty) . and h is a harmonic field, i. e. 

(3.4) div v = 0, div h = 0, V/i = 0. 

In particular the subspace of all harmonic tensor fields in L 2 (T£M) is of 
finite dimensional: 

(3.5) H{T r s M) = {h G L 2 (T s r M)| V/i = 0, div ft = 0}, dim H < oo. 

Remark 3.1. The above orthogonal decomposition theorem implies that L 2 (E) 
(E = TgM) can be decomposed into 

L 2 (£) =(?(£■)© L 2 D (E) for 0M = 0, 

L 2 {E) = G(E) © H{E) © for ill compact. 

Here H is as in (|3.5j) . and 

G(£) = {^i 2 (£) | i; = V^(^ei/ 1 (r;_ 1 M)}, 

L 2 D {E) L 2 (£T) | div v = 0}, 

4(£) = {,e4(B) I V^0}. 

They are orthogonal to each other: 

L 2 D (E)±G(E), L 2 N (E)+H(E), G(E)±H(E). 

Remark 3.2. The dimension of the harmonic space H(E) is related with the 
bundle structure of E — TJM. It is conjectured that 

dimi? = k = the degree of freedom of E. 

Namely k is the integer that E can be decomposed into the Whitney sum of a 
fc-dimensional trivial bundle E k = Mxl' and a nontrivial bundle E%, i.e. 

E = Ei © E k . 

Proof of Theorem \3.1\ We proceed in several steps as follows. 

Step 1 Proof of Assertion (1). Let u G L 2 (E), E = T r s M (r + s > 1). 
Consider the equation 

(3.6) Ap — div u in M, 
where A is the Laplace operator defined by 

(3.7) A = divV. 
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Without loss of generality, we only consider the case where div u G E = TJ _1 M. 
It is clear that if the equation (|3.6p has a solution ip 6 if 1 (£7), then by (|3.7[) . the 
following vector field must be divergence-free 

(3.8) v = u-Vip e L 2 (E). 
Moreover, by (|3.1|) we have 

(3.9) (v,V<p)=0. 

Namely v and \7<p are orthogonal. Therefore, the orthogonal decomposition u = 
v + Vip follows from ((378]) and (|3J)l . 

It suffices then to prove that Q3.6j) has a weak solution 6 H 1 {E): 

(3.10) (V^ - w, V-0) = y^eH\E). 
To this end, let 

H = H\E)\H, 
H = {ip e H 1 ^)] V^> = 0}. 
Then we define a linear operator G : H — > H by 

(3.11) {G<p,il>) = (Vy, Vt/O V-0 G ff. 

It is clear that the linear operator G : H — > H is bounded, weakly continuous, and 

(3.12) (G^) = (V^,V^) = |M| 2 . 
Based on Theorem [221 f° r an y / G H , the equation 

= / in M 



has a weak solution ip £ H. Hence for / = divit the equation (|3 .6[) has a solution, 
and Assertion (1) is proved. In fact the solution of (|3.6[) is unique. We remark that 
by the Poincare inequality, for the space H = H 1 (E) \ {t/j\Wijj = 0}, (|3.12j) is an 
equivalent norm of H . In addition, by Theorem 12.11 the weak formulation (|3.10|) 
for (|3.6p is well-defined. 

Step 2 Proof of Assertion (2). Based on Assertion (1), we have 

H k (E) = H k D ®G k , 
L 2 {E) = L 2 D ®G, 

where 

H k ) = {ueH k (E) | divu = 0}, 
G k = {ue H k (E) | u^Wtp}. 

Define an operator A : Hp (E) — > L 2 D (E) by 

(3.13) Au = PAu, 

where P : L 2 (E) — > L 2 D (E) is the canonical orthogonal projection. 
We know that the Laplace operator A can be expressed as 

(3.14) A = divV = ^lW^ + * ) 
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where B is the lower order derivative linear operator. Since M is compact, the 
Sobolev embeddings H 2 (E) °-> H 1 (E) L 2 (E) are compact, which implies that 
the lower order derivative operator 

B : H 2 (M,R N ) -> L 2 (M,R N ) is compact, 

where the integer N is the dimension of the tensor bundle E. According to the 
elliptic operator theory, the elliptic operator in (|3.14[) 

A = '/ l -j^ iJ? ■■ H 2 (M,R N ) -+ L 2 (M,R N ) 

is a linear homeomorphism. Therefore the operator in (|3.14p is a linear completely 
continuous field 

A : H 2 (E) L 2 (E), 
which implies that the operator of (|3.13[) is also a linear completely continuous field: 

A = PA: H 2 D {E) -> L 2 D (E). 

By the the spectral theorem of completely continuous fields [51 15] , the space 

H = {u G H 2 D {E)\ Au = 0} 

is finite dimensional, and is the eigenspace of the eigenvalue A = 0. By Theorem 
EU for u G H 



(Au,u)^—gdx = / (Au, u)yj—gdx (by divu = 0) 

' m Ju 

= — I (V«, Wu)yf—gdx 

JAI 

=0 (by Am = 0). 

It follows that 

!ieifoV« = H = H, 
where H is the harmonic space as in (|3.5p . Thus we have 

L 2 D (E) = H(BL%(E), 

L 2 N (E) = {u G L 2 D (E)\ Vu^O}. 
The proof of Theorem 13. II is complete. □ 

3.2. Uniqueness of the orthogonal decomposition. In Theorem 13. 11 a tensor 
field u E L 2 (T^M) with r + s > 1 can be orthogonally decomposed into 

u = Vtp + v for dM = 0, 

^ u = V</5 + v + h for M compact. 

Now we address the uniqueness problem of the decomposition (|3.15l) . In fact, if 
u is a vector field or a co-vector field, i.e. 

u G L 2 (TM) or u G L 2 (T*M), 

then the decomposition of (|3.15|) is unique. 

We can see that if u G L 2 (TJM) with r + s > 2, then there are different types of 
the decompositions of (|3.15[) . For example, for u G ^(XfM), the local expression 
of it is given by 

U = {Uij(x)}. 
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In this case, u has two types of decompositions: 

(3.16) Uij = D,(pj + Vij, D l Vij = 0, 

(3.17) u^ = Djipi + Wij, D l Wij = 0. 

It is easy to see that if ^ Uji then both (|3.16[) and (|3.17l) are two different 
decompositions of Uij . Namely 

{%} 7^ { w ij}> (^1) ' ' ' )<Pn) (i>l, ■ ■ ■ >i>n)- 

If Uij = Uji is symmetric, u can be orthogonally decomposed into the following 
two forms: 

= + Di<pj, D l Vij = 0, 
= Wij + Djtpi, D l Wij = 0, 

and <p and ip satisfy 

(3.18) Aipj = D k u kj , 



k ■ 



(3.19) Aipj = D K Uj 

By u^ — Uji we have D k Ukj = D k Uj k - Hence, (13.181) and (|3.19p are the same, and 
ip = ip. Therefore, the symmetric tensors Uij can be written as 

(3.20) u^ = + Diipj, D l v l3 = 0, 

(3.21) u^ = Wij + Dj(fi, D j Wij = 0. 

From (|3.20[) - (|3.21[1 we can deduce the following theorem. 

Theorem 3.2. Let u G L^T^M) be symmetric, i.e. d the first Betti 

number f3\(M) = for M . Then the following assertions hold true: 

(1) u has a unique orthogonal decomposition if and only if there is a scalar 
function (p € H 2 (M) such that u can be expressed as 

(3 22) = + DiDjtp, 

= Vji, D l Vij = 0. 

(2) u can be orthogonally decomposed in the form of \3. 22]) if and only if u^ 
satisfy 

where R k — g kl Rij and Rij are the Ricci curvature tensors. 

(3) If Vij in $3.20\) is symmetric: Vij — Vji, then u can be expressed by (3.22\) . 

Proof. We only need to prove Assertions (2) and (3). 
We first prove Assertion (2). It follows from (13.20[) that 

,3.24, A (DVi) Bp 0A Vl 



dxi dx % dx 3 dx 

where A = D k D k . By the Weitzenbock formula [3], 

(3.25) A<pi = -(Sd + d% 4 - R k ip k , 
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and (Sd + dS) is the Laplace-Beltrami operator. We know that for ui = ipidx 1 

dip 



dui = tfi 

dSco = V(Ayj) tfi 
where V is the gradient operator, and 



dx l ' 
dip 

dx % ' 



1 d ( i—. a d( P 



Namely 



d dip 
(Sd + dS)cpi = —Kip Vi = ^ 



Hence, we infer from (|3.25p that 

Thus, by (pT24|) and (|3T26|) . we obtain that 

holds true if and only if the tensor -0 = (<^i, • • • , ip n ) in (I3.20|) is a gradient ip = Vip. 
Assertion (2) is proven. 

Now we verify Assertion (3). Since Vij in (13.20[) is symmetric, we have 

(3.27) Di<pj = Djcpi. 

Note that 



where is the Levi-Civita connection, and — Tj t . We infer then from (|3.27p 
that 

^ ^ dx i dxi 

By assumption, the 1-dimensional homology of M is zero, 

Hi (AO =0, 

and it follows from the de Rham theorem and (|3.29[) that 

dip 

fk = q-T) 
OX K 

for some scalar function ip. Thus Assertion (3) follows and the proof is complete. □ 



Remark 3.3. The conclusions of Theorem 13.21 are also valid for second-order 
contra- variant symmetric tensors u — {u lJ }, and the decomposition is given as 
follows: 

u ij =v ij +g ik g jl D k D l ip, 

D iV 13 = 0, v ij = v j \ ip £ H 2 (M). 
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4. Variational Principle for Functionals of Riemannian Metric 

4.1. General theory. Hereafter we always assume that M is a closed manifold. 
A Riemannian metric G on M is a mapping 

(4.1) G: M -+T$M = T*M®T*M, 

which is symmetric and nondegenerate, i.e., in a local coordinate (v, x),G can be 
expressed as 

(4.2) G = {g ij [x)} with </,, '/„. 
and the matrix (g%j) is invertible on M: 

(4-3) (gV) = ( gij )-\ 

If we regard a Riemannian metric G = {gij} as a tensor field on manifold M, 
then the set of all metrics G = {gij} on M constitutes a topological space, called 
the space of Riemannian metrics on M. We denote 

(4.4) G" 1 = {g ij } : M -> T 2 M = TM <g) TM. 
The space for Riemannian metrics on M is denied by 

(4.5) W m ' 2 (M,g) = {g I G € W m ' 2 {T%M),G- 1 e W m ' 2 (T 2 M), 

G is the Riemannian metric on M as in (|4.2|) }. 

which is a metric space, but not a Banach space. However, it is a subspace of the 
direct sum of two Sobolev spaces W m ' 2 {T$M) and W m < 2 {T$M): 

W m ' 2 (M,g) C W m ' 2 {T°M) 8 W m 2 {T 2 M). 

A functional defined on W m ' 2 (M, g) : 

(4.6) F : W m ' 2 (M,g) -> K 

is called the functional of Riemannian metric. Usually, the functional (|4.6p can be 
expressed as 

(4.7) F( 9ij )= f f(x, 9ij ,--- ,d m giJ )^—gdx. 

J M 

Since {g^) is the inverse of (gij), we have 

(4.8) gij = — x the cofactor of g y . 

Therefore, F(gij) in (|4.7p also depends on g l ° , i.e. putting (|4.8I) in (|4.7|) we get 

(4.9) F(g i i)= [ /;,■.</"...../)-,/, ; •,—/,-. 

We note that although Ty m,2 (M,<7) is not a linear space, for a given element 
<7ij G W n,p (M, and any symmetric tensor Xij and there is a number Ao > 
such that 

i S« + XXa € V^ m ^ 2 (M, 5 ) VO < |A| < Ao, 

' 4 " ' 1 ' ' g ij + XX 11 e W m - 2 (M, g) VO < |A| < A . 
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Due to (|4.10[) . we can define the derivative operators of the functional F, which 
are also called the Euler-Lagrange operators of F, as follows 



(4.11) 6*F : W m ' 2 (M,g)^W~ m - 2 (T 2 M), 

(4.12) S*F : W m - 2 {M,g)^W- m - 2 {T°M), 

where W~ m ' 2 {E) is the dual space of W m ' 2 (E), and S*F,5*F are given by 

(4-13) (6*F( 9ij ),X) = ^F{ 9ij + XX l3 )\ x=0 , 

(4.14) (5*F(gV),X) = ^-F(g^ + AA^)| A=0 . 



For any given metric gij £ W m,2 (M, g), the value of 8*F and S*F at gij are 
second-order contra-variant and covariant tensor fields respectively, i.e. 

5*F( 9ij ) : M -t TM ® TM, 
5*F( 9ij ) :M^T*M® T*M. 

Moreover, the equations 

(4.15) S*F( 9ij ) = 0, 

(4.16) 8*F( 9ij ) = 0, 

are called the Euler-Lagrange equations of F, and the solutions of (|4. 1 5[) and (|4.16[) 
are called the extremum points or critical points of F. 

Theorem 4.1. Let F be the junctionals defined by |^.6'[ ) and Then the 

following assertions hold true: 

(1) For any g^ e W m ' (M,g),5*F( 9 ij) and S*F(gij) are symmetric tensor 
fields. 

(2) If {g^} £ W rn ' 2 (M, g) is the extremum point of F, then {g^} is also an 
extremum point of F, i.e. {g^} satisfies . _? 5| ) and J^._?6'| ) if and only if 
{g l °} satisfies h4.15ty and 14.16\ ). 

(3) 6*F and S*F have the following relation 

(6*F( 9ij )) kl = -g kr g ls (6*F(g l3 )) rs , 
where (5* F) kl and (S*F)ki are the components of 5* F and 5*F respectively. 
Proof. We only need to verify Assertion (3). Noting that 

9ik9 kj =S{, 

we have the variational relation 

K9ik9 kj ) = 9tkSg k3 + g k3 5gik = 0. 

It implies that 

(4.17) 6g kl = - 9 ki 9 lj 6 9ij . 
In addition, in flUTS} and (|414|) . 

XXij = Sg %3 , XX 11 = Sg 13 , X^O small. 
Therefore, by (14. 17)) we get 

((5*F) kl ,5g kl ) = -{{5,F) kU g kl g l ^ 9lJ ) = (-g ki g lj (6*F) kl ,S 9ij ) = ((5*F)»,S 9ij ). 
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Hence 

(5*F)« =-g ki g lj (6*F) kl . 
Thus Assertion (3) follows and the proof is complete. □ 

4.2. Scalar potential theorem for constraint variations. We know that the 
critical points of the functional F in (|4.6[) are the solution 

(4.18) SF( 9ij ) = 0, 
in the following sense 

(4.19) (SF( 9ij ),X) =-^F(9 ij + AX«)| A=0 

= / {5F{ gij ))uX kl y/=g~dx 

JM 

=0 yX kl = X lk in L 2 {E), 

where E = TM ® TM . Hence, the critical points of functionals of Riemannian 
metrics are not solutions of (|4.18[) in the usual sense. 

It is easy to see that L (TM ® TM) can be orthogonally decomposed into the 
direct sum of the symmetric and contra-symmetric spaces, i.e. 

L 2 (E)=L 2 S (E)®L 2 C (E), 

(4.20) L 2 (E) = {ue L 2 (E)\ Uij = Uji }, 

L 2 c (E) = {veL 2 (E)\ v lJ =-v lJ }. 

Since 8F is symmetric, by (|4.20j) the extremum points {gij} of F satisfy more 
general equality 

(4.21) (5F( 9ij ),X) = VX = {X tJ } e L 2 (E). 

Thus, we can say that the extremum points of functionals of the Riemannian 
metrics are solutions of (|4.18j) in the usual sense of (|4.21|) , or are zero points of the 
variational operators 

SF: W m > 2 (M,g) ^W~ m > 2 (E). 

Now we consider the variations of F under the divergence-free constraint. In this 
case, the Euler-Lagrangian equations with symmetric divergence-free constraints 
are equivalent to the Euler-Lagrangian equations with general divergence-free con- 
straints. Hence we have the following definition. 

Definition 4.1. Let F : W m ' 2 (M, g) — > M be afunctional of Riemannian metric. A 
metric tensor {g^} £ W m ' 2 (M, g) is called an extremum point of F with divergence- 
free constraint, if {gij} satisfies 

(4.22) (5F( 9ij ),X)=0 VX = {X ij }cL 2 D (E), 
where L 2 D (E) is the space of all divergence-free tensors: 

L 2 D (E) ={Ie L 2 {E) | div X = 0}. 

It is clear that an extremum point satisfying (|4.22p is not a solution of (|4.18l) . 
Instead, we have the scalar potential theorem for the extremum points of divergence 
free constraint (|4.22|) . which is based on the orthogonal decomposition theorems. 
This result is also crucial for the gravitational field equations and the theory of 
dark matter and dark energy developed later. 
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Theorem 4.2 (Scalar Potential Theorem). Assume that the first Betti number of 
M is zero, i.e. f3\{M) = 0. Let F be a functional of the Riemannian metric. Then 
there is a <p e H 2 (M) such that the extremum points {gij} of F with divergence-free 
constraint satisfy 

(4.23) (SF( 9ij )) k i = D k Dnp. 

Proof. Let {gtj} be an extremum point of F under the constraint (|4.22|) . Namely, 
SF(gij) satisfies 

rfci 



(4.24) (5F( 9ij )) kl X kl ^dx = VX = {X kl } with D k X kl = 0. 

J M 

By Theorem 13. 1[ S*F(gij) can be orthogonally decomposed as 

(4.25) {5F{g l0 )) u = v kl + D k ^, D k v kl = 0. 
By Theorem EU for any D k X kl = 0, 

kl\ / 7~i vkl / ~ 7_ / l r} vkl 



(4.26) (D k ^,X M )= D k iP l X kl ^g~dx= - ^D k X kl ^jdx = 0. 

Jm Jm 

Therefore it follows from (|4T24>(|4T2l))) that 

'kl I -J™ _ n wr>. vkl _ n 



(4.27) / v kl X kl ^dx = VD k X k 

Jm 

Let X kl = g kl g l] v rj . Since 

D kgij = D k g ij = 0, 

we have 

D k X kl = />,;//'//'r, ; i = g l ^g ik D kVij ) = g' 1 !)',,, = 
thanks to D l v l3 = 0. Inserting X kl = g ki g l] v t] into (|4.27|) leads to 

2 _ / n ki„lj, 



M\l2= g kl g t3 v k iv ijv ^gdx = 0, 

J M 

which implies that v = 0. Thus, (|4.25|) becomes 
(4.28) (SF( 9ij )) kl =D k tl>t. 

By Theorem 14. 11 Ji* 1 is symmetric. Hence we have 

D k ij)i = Di^p k . 

It follows from ([3~2"gjl that 

By assumption, the first Betti number of M is zero, i.e. the 1-dimensional homology 
of M is zero: H\(M) = 0. It follows from the de Rham theorem that if 

d (^te*)=Qj*-0) dx'Adi fc = 0, 

then there exists a scalar function <£> such that 

^ = ^rrdx k — ip k dx k . 
ox" 

Thus, we infer from (|4.29l) that 

= for some V G H 2 (M). 
ox 1 
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Therefore we get (|4.23p from (|4.28[) . The theorem is proved. □ 

If the first Betti number f3\{M) ^ 0, then there are N — f5\{M) number of 
1-forms: 

(4.30) ujj — ip^ k dx k £ H^(M) for 1 < j < N, 

which constitute a basis of the 1-dimensional de Rham homology H\(M). We know 
that the components of uij are co- vector fields: 

(4.31) ip j = ■ ■ ■ , ijp n ) £ H\T*M) for 1 < j < N, 
which possess the following properties: 

dx l dx k — J — i 

or equivalently, 

Dii>{ = D kM for 1 < 3 < N- 

Namely, Vip 11 £ L 2 (T*M ® T*M) are symmetric second-order contra-variant ten- 
sors. Hence Theorem 14.21 can be extended to the non- vanishing first Betti number 
case as follows. 

Theorem 4.3. Let the first Betti number pi(M) ^ for M. Then for the func- 
tional F of Riemannian metrics, the extremum points {gij} of F with the constraint 
\4 ■ 22[) satisfy the equations 

N 

(4.32) (5F( gij )) kl = D k Dicp + J2ajDki4> 

3=1 

where N = (3i(M), ctj are constants, ip £ H 2 {M), and the tensors ^ = • • • , ip^) 
H 1 (T*M) are as given by f£3?P . 

The proof of Theorem 14.31 is similar to Theorem 14.21 and is omitted here. 

Remark 4.1. By the Hodge decomposition theory, the 1-forms ujj in (|4. 30(1 are 
harmonic: 

duij = 0, Sojj = for l<j<N, 
which implies that the tensors ipi in (|4.32[) satisfy 

(4.33) (Sd + dS)^ j = for 1 < j < N. 
According to the Weitzenbock formula (|3.25|) . we obtain from (|4. 3311 that 

(4.34) F) k D k ^\ = -Rftpi for 1 < j < N, 
for ^ = in (14321 . 



Remark 4.2. Theorem [42] is derived for deriving new gravitational field equation 
in the next section for explaining the phenomena of dark matter and dark energy. 
The condition that j3i (M) = means that any loops in the manifold M can shrink 
to a point. Obviously, our universe can be considered as a 4-dimensional manifold 
satisfying this condition. 
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Part 2. Physics 

5. Gravitational Field Equations 

5.1. Einstein-Hilbert functional. The general theory of relativity is based on 
three basic principles: the principle of equivalence, the principle of general relativ- 
ity, and the principle of Lagrangian dynamics. The first two principles tell us that 
the spatial and temporal world is a 4-dimensional Riemannian manifold (M, gij), 
where the metric {gij} represents gravitational potential, and the third principle de- 
termines that the Riemannian metric {g^ } is an extremum point of the Lagrangian 
action, which is the Einstein-Hilbert functional. 

Let (M.gij) be an rt-dimensional Riemannian manifold. The Einstein-Hilbert 
functional 

(5.1) F : W 2 ' 2 (M, g) -> R 
is defined by 

(5.2) F( 9ij ) = J (r + ^V%) V^9dx, 

where W 2 ' 2 {M,g) is defined by (|I3]> , R = g kl R kl and R kl are the scalar and the 
Ricci curvatures, S%j is the stress tensor, G is the gravitational constant, and c is 
the speed of light. 

The Euler-Lagrangian of the Einstein-Hilbert functional F is given by 

(5-3) 5F( 9ij ) = Rij - ± 9ij R + 

where Ty is the energy-momentum tensor given by 

Tij = Sij-\ 9ij S + g kl ^, S = g kl S kl , 
and the Ricci curvature tensor Rij is given by 

(5.4) RlJ = Vf^+ ^ m <r2,hl ' ,2,IL > 



13 2 \dx l dx : > dx k dx l dx^dx k dx l dx l 

By (|5.3p - (|5.5p . the Euler-Lagrangian SF(gij) of the Einstein-Hilbert functional is a 
second order differential operator on {gij}, and SF(gij) is symmetric. 

5.2. Einstein field equations. The General Theory of Relativity consists of two 
main conclusions: 

1) The space-time of our world is a 4-dimensional Riemannian manifold (M 4 , gij), 
and the metric {gij} represents gravitational potential. 

2) The metric {g^} is the extremum point of the Einstein-Hilbert functional 
(|5.2p . In other words, gravitational field theory obeys the principle of La- 
grange dynamics. 



^The matter tensor is included here as well. 
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(5.8) Rij - -gijR — T ij 



The principle of Lagrange dynamics is a universal principle, stated as: 

Principle of Lagrange Dynamics. For any physical system, there are a set of 
state functions 

u = ■ • • ,u N ), 

which describe the state of this system, and there exists a functional L of u, called 
the Lagrange action: 

(5.6) L{u)= [ [ £(u,Du,--- ,D m u)dxdt, 

Jo Jn 

such that the state u is an extremum point of L. Usually the function £ in H5. 6}) is 
called the Lagrangian density. 

Based on this principle, the gravitational held equations are the Euler-Lagrange 
equations of the Einstein-Hilbert functional: 

(5.7) 6F( 9ij ) = 0, 
which are the classical Einstein held equations: 

1 „ 8ttG, 
T 

By the Bianchi identities, the left hand side of (|5.8[) is divergence-free, i.e. 

(5.9) D\R l3 - ± 9ij R) = 0. 

Therefore it is required in the general theory of relativity that the energy-momentum 
tensor {Tij} in (|5.8|) satisfies the following energy-momentum conservation law: 

(5.10) D l T %j = for 1 < j < n. 

5.3. New gravitational field equations. Motivated by the mystery of dark en- 
ergy and dark matter and the other difficulties encountered by the Einstein held 
equations as mentioned in Introduction, we introduce in this section a new set of 
held equations, still obeying the three basic principles of the General Theory of 
Relativity. 

Our key observation is a well-known fact that the Riemannian metric gij is 
divergence-free. This suggests us two important postulates for deriving the new 
gravitational held equations: 

• The energy-momentum tensor of matter need not to be divergence-free due 
to the presence of dark energy and dark matter; and 

• The field equation obeys the Euler-Lagrange equation of the Einstein- 
Hilbert functional under the natural divergence- free constraint. 

Under these two postulates, by the Scalar Potential Theorem, Theorem 14.21 if the 
Riemannian metric {gij} is an extremum point of the Einstein-Hilbert functional 
(|5.2p with the divergence- free constraint (|4.22p . then the gravitational field equa- 
tions are taken in the following form: 

1 &7rG 

(5.11) R^ - -g i:j R = — —Tij - D l D j tp, 

where 99 € H 2 (M) is called the scalar potential. We infer from (|5.9I) that the 
conservation laws for (|5.11|) arc as follows 

(5.12) div (DiDjip + ^fTij) = 0. 
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Using the contraction with g 1 ^ in (|5.11|) , we have 

(5.13) R=^.T + 9, 

where 

/• .'/"/;.,• n"ihi),^. 

represent respectively the energy-momentum density and the scalar potential den- 
sity. Physically this scalar potential density $ represents potential energy caused 
by the non-uniform distribution of matter in the universe. One important property 
of this scalar potential is 



(5.14) / <Py/^)dx = 0, 

Jm 

which is due to the integration by parts formula in Theorem 12.11 This formula 
demonstrates clearly that the negative part of this quantity $ represents the dark 
matter, which produces attraction, and the positive part represents the dark energy, 
which drives the acceleration of expanding galaxies. We shall address this important 
issue in the next section. 

5.4. Field equations for closed universe. The topological structure of closed 
universe is given by 

(5.15) M = S 1 x S 3 , 

where S 1 is the time circle and S 3 is the 3-dimensional sphere representing the 
space. We note that the radius R of S 3 depends on time t G S 1 , 

R = R(t), teS 1 , 

and the minimum time to, 

to = min R(t) 

is the initial time of the Big Bang. 

For a closed universe as (|5.15|) . by Theorem 14. 2 [ the gravitational field equations 
are in the form 

1 &ttG 

(5 - 16) A^+ 5 lfc i?^ fc =0, 
Ditpj = Djijii, 

where A = D k Dk,(fi the scalar potential, ip = (tpo, ^1,^2,^3) the vector potential, 
and a is a constant. The conservation laws of (|5 . 16[) are as follows 

1 A (d<p\ , 8^G ofcr 



(5-17) A^ i = -A^j + - 3 -^T fc ,. 

6. Interaction in A Central Gravitational Field 

6.1. Schwarzschild solution. We know that the metric of a central gravitational 
field is in a diagonal form pQ: 

(6.1) ds 2 = g m c 2 dt 2 + 9ll dr 2 + r 2 (d9 2 + sin 2 6dip 2 ), 
and physically goo is given by 

(6.2) 5 oo = - U + ~^l> 
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where ip is the Newton gravitational potential; see among others [I] . 

If the central matter field has total mass M and radius ro, then for r > ro, the 
metric (|6.1I) is the well known Schwarzschild solution for the Einstein field equations 
(|5.8p . and is given by 

(6.3) ds 2 = - (l - c 2 dt 2 + dr ' MG + r 2 (d8 2 + sin 2 9d<p 2 ). 
We derive from (|6.2|) and (|6.3j) the classical Newton gravitational potential 

(6.4) V = - — • 

r 

6.2. New gravitational interaction model. We now consider the metric de- 
termined by the new field equations (|5.11[) . from which we derive a gravitational 
potential formula replacing (|6.4p . 

Equations (|5.11|) can be equivalently expressed as 

-gijT) - (DiDjip - -i 
where 

T = g kl T kh <P = g kl D k D Wl 

For the central matter field with total mass M and radius ro , by the Schwarzschild 
assumption, for r > ro, there exists no matter, i.e. 

(6.6) Tij = 0. 

Therefore the conservation laws of 

dip 
dx 1 



(6-5) Rij = {T l0 - - 9ij T) - (DiDj <p - - 9lJ H) , 



(6.7) A(^)=0 forfc = 0,l,2,3. 



The tensors gij in (|6.ip are written as 

,„ oN .9oo = -e n , 5n = e°, g 22 = r 2 , g 33 = r 2 sin 2 9, 

u = u{r), v = v(r). 

Noting that the central field is spherically symmetric, we assume that 

(6.9) tp — <p(r) is independent of t, 9, ip. 

2MG 

(6.10) r >__. 



For the metric (|6.7|) . the non-zero components of the Levi-Civita connection are 
as follows 

1 00 



(6.11) I* 



33 



r 2 

1 33 



= -e u - v u' 
2 


V 1 - 
1 n 




r 1 - 

1 22 — 


—re 


= -re~ v sin 2 9, 


1 10 — 


2 ' 


F 2 - 

1 12 — 


1 

r ' 


= — sin 9 cos 9, 


p3 _ 
1 13 — 


1 

* 

r 


r 3 - 

1 23 — 


cos 
sin 9 



Hence the Ricci tensor 



D _ tfc _ 2 J _|_ pfc pr pre pr 

3 ~ dxi dx k jk ij kr 
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are given by 



R 



no 



u" v! v! , . „ 



(6.12) 



„ u v u ' , ,s 
R,i = --- + 1 (u>-v>), 



1-e* + -(«'-!/) 



i?22 = E 

R33 = -R22 sin^ 9, 
Rj =0 Vi ^ ./• 
Furthermore, we infer from (16.71) . (|6.9[) and (|6.1ip that 

1 



(6.13) 





2*** 


D D a ip - 


1 


D\Dxip - 




D 2 D 2 (fi - 


^322* 


D 3 D 3 (p - 


^933^ 



1 4 
" iii 1 \ 1 
if - -(u +v - -)(p 



<p"-\w + u' + -W 

2 r 



ip" + ~(u' - v')ip 



1 



Thus, by (|6.12[) and (|6. 13[) . the equations (|6.5[) are as follows 
(6.14) 
(6.15) 
(6.16) 



u" + — + %{u' v>) = y" - Uu' +v'- -)tp>, 
r 2 2 r 

u" - — + %{u'- v') = -v" + \{v! + v> + V, 



u' - t/ + -(1 - e") = r(<p" + l(u' - v')cp'). 
r '2 



6.3. Consistency. We need to consider the existence and uniqueness of solutions 
of the equations (|6.14[) - (|6.16|) . First, in the vacuum case, the classical Einstein 
equations are in the form 



(6.17) 



R kk =0 for k = 1,2,3, 



two of which are independent. The system contains two unknown functions, and 
therefore for a given initial value (as u' is basic in (|6.17p 'l: 



(6.18) 



u'(r ) = (Ti, u(r ) = cr 2 , r o > 0, 



the problem (|6.17[) with (|6.18[) has a unique solution, which is the Schwarzschild 
solution 



S-2 



1 



In 1 



£1 



where 



ei =r (l-e- CT2 ), e 2 = 4 1- 



£1 
ro 



(Tl. 
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Now if wc consider the influence of the cosmic microwave background (CMB) 
for the central fields, then we should add a constant energy density in equations 

f-g 00 p 0\ 





\ 0/ 



(6.19) 



(Tij) 



Namely, 
(6.20) 



_ 4ttG 

-KOO — — 7-g00P, 



Rn — 0, R22 — 0, 



where p is the density of the microwave background, whose value is p = 4 x 
10~ 31 kg/m 3 . Then it is readily to see that the problem (|6.20j) with (|6.18j) has 
no solution. In fact, the divergence-free equation D l Tij — yields that 

r? r o = \u'p = 0, 



which implies that v! = 0. Hence R 0Q — 0, a contradiction to (|6.20p . Furthermore, 
if we regard p as an unknown function, then the equations (|6.20p still have no 
solutions. 

On the other hand, the new gravitational field equations (|6.14p - (l6.16[) are solvable 
for the microwave background as the number of unknowns are the same as the 
number of independent equations. 

Equations (16.14j) - (|6.16|) have the following equivalent form: 



(6.21) 



1 u' 
r + 2 
1 , ., 



1) 



r 

1 / , 1 , 
-ipv + -u 

2 r 



i(^-l)-^V, 



equipped with the following initial values: 

(6.22) u'(r ) = a x , v(r ) = a 2 , <p'(r ) = a 3 , r Q > 0. 

It is classical that (|6.21j) with (|6.22j) possesses a unique local solution. In fact, we 
can prove that the solution exists for all r > tq. 



6.4. Gravitational interaction. We now derive the gravitational interaction for- 
mula from the basic model (|6.14[) - (|6.16|) . 
First, we infer from (f6TT4|) - (f6TT6|) that 



(6.23) 

Consequently, 
(6.24) 



u + v = 



rip 



r r ' 



-(e v -l) + r<p" 



r 1 

2<P r 



1) 



rip 
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By ((621) and <|6?Tj) . we have 
(6.25) 

As the interaction force F is given by 



V;=y(^-1). 



it follows from ([6723]) and (|6~25|) that 

™„2 



(6.26) 



3.27) 



F 



1 



2 r 



rp 



Of course, the following energy balance and conservation law hold true as well: 

iirG r 



(6.28) 



R 



3 «+V$dr = 0. 



where i? is the scalar curvature and $ = g kl DkDiip. Equation (|6.27|) is derived by 
solving ip" using (|6.28|) . Namely, for r > ro where T = 0, 



i? = $ = g M D k Di<p 



2 r 



6.5. Simplified formulas. We now consider the region: ro < r < r\. Physically, 
we have 

(6.29) |^|,|^"| «1. 

Hence u and v in (|6.26l) can be replaced by the Schwarzschild solution: 



(6.30) 



uq = In ( 1 — 



«o = - In ( 1 - 



5 = 



2GM 



As S/r is small for r large, by (|6.29j) . the formula (|6.26j) can be expressed as 



~2 r 



(6.31) 



F = 



'12 ~ ^ r 



This is the interactive force in a central symmetric held. The first term in the 
parenthesis is the Newton gravity term, and the added second term —rip" is the 
scalar potential energy density, representing the dark matter and dark energy. 

In addition, replacing u and v in (|6.27j) by the Schwarzschild solution (|6.30j) . we 
derive the following approximate formula: 

'2 + 5 



(6.32) 



9 



ip 1 - R. 



Consequently we infer from (|6.31j) that 



(6.33) 



F = mMG 



1 1 



--o-7 2 + - W + — 



Ri 



R = $ for r > r . 



The first term is the classical Newton gravitation, the second term is the coupling 
interaction between matter and the scalar potential <p, and the third term is the 
interaction generated by the scalar potential energy density $. In this formula, the 
negative and positive values of each term represent respectively the attracting and 
repelling forces. 
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Integrating (|6.32p yields (omitting e s ' r ) 



3.34) 



<P 



-e 2 r 



l Rdr, 



where £2 is a free parameter. Hence the interaction force F is approximated by 



(6.35) 



F = mMG 



1 



Rr 1 



— - + 2 + - + — + - 2 + - 



l Rdr 



where e = £2<5~ 1 , R = for r > ro, and 5 = 2MG/c 2 . We note that based on 
(16.281) . for r > r , R is balanced by and the conservation of $ suggests that i? 
behaviors like r~ 2 as r sufficiently large. Hence for r large, the second term in the 
right hand side of f|6 . 35[) must be dominate and positive, indicating the existence 
of dark energy. 

We note that the scalar curvature is infinite at r = 0: R(0) = 00. Also R contains 
two free parameters determined by u' and v respectively. Hence if we take a first 
order approximation as 

(6.36) R = -£1 + — for r < r < n = 10 21 fcm, 

r 

where e\ and £0 are free yet to be determined parameters. Then we deduce from 
(f6T34l) and (j6~36| that 



£"0 



Therefore, 
(6.37) 



F = mMG 



I c 2 U 

99 = — £or — £ir + — 



-\-^--+£l + (eS + £l 5- l )r + 2er 2 
v l 2 r 



Physically it is natural to choose 

£0 > 0, £1 > 0, £ > 0. 
Also, £1 and 2£2^ _1 r 2 are much smaller than (e5 + £i^ _1 )r for r < r\. Hence 

1 fco 



(6.38) 



F = mMG 



where fco and k\ can be estimated using the Rubin law of rotating galaxy and the 
acceleration of the expanding galaxies: 



(6.39) 



fcn = 4 x 10" 



fci = icr a 'fcm 



,-3 



We emphasize here that the formula (|6.38[) is only a simple approximation for 
illustrating some features of both dark matter and dark energy. 



7. Theory of Dark Matter and Dark Energy 

7.1. Dark matter and dark energy. Dark matter and dark energy are two of 
most remarkable discoveries in astronomy in recent years, and they are introduced 
to explain the acceleration of the expanding galaxies. In spite of many attempts 
and theories, the mystery remains. As mentioned earlier, this article is an attempt 
to develop a unified theory for the dark matter and dark energy. 

A strong support to the existence of dark matter is the Rubin law for galactic 
rotational velocity, which amounts to saying that most stars in spiral galaxies orbit 
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at roughly the same speed. Namely, the orbital velocity v(r) of the stars located at 
radius r from the center of galaxies is almost a constant: 

(7.1) v(r) = constant for a given galaxy. 

Typical galactic rotation curves [3] are illustrated by Figure ITTTt a). where the verti- 
cal axis represents the velocity (km/s), and the horizontal axis is the distance from 
the galaxy center (extending to the galaxy radius). 
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(a) 



(b) 



Figure 7.1. (a) Typical galactic rotation curve by Rubin, and (b) 
theoretic curve by the Newton gravitation law. 



However, observational evidence shows discrepancies between the mass of large 
astronomical objects determined from their gravitational effects, and the mass cal- 
culated from the visible matter they contain, and Figure 17.11 (b) gives a calculated 
curve. The missing mass suggests the presence of dark matter in the universe. 

In astronomy and cosmology, dark energy is a hypothetical form of energy, which 
spherically symmetrically permeates all of space and tends to accelerate the expan- 
sion of the galaxies. 

The High-Z Supernova Search Team in 1998 and the Supernova Cosmology 
Project in 1998 published their observations which reveal that the expansion of 
the galaxies is accelerating. In 2011, a survey of more than 2 x 10 5 galaxies from 
Austrian astronomers confirmed the fact. Thus, the existence of dark energy is 
accepted by most astrophysicist. 

7.2. Nature of dark matter and dark energy. With the new gravitational 
field equation with the scalar potential energy, and we are now in position to derive 
the nature of the dark matter and dark energy. More precisely, using the revised 
Newton formula derived from the new field equations: 



(7.2) 



F = m 'M G ( -i - — + k 1 r 



we determine an approximation of the constants fco,fci, based on the Rubin law 
and the acceleration of expanding galaxies. 

First, let M r be the total mass in the ball with radius r of the galaxy, and V 
be the constant galactic rotation velocity. By the force equilibrium, we infer from 
CL2]) that 
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which implies that 

^ M - = ^ l + fco r- fcl r3 - 

This matches the observed mass distribution formula of the galaxy, which can 
explain the Rubin law (|7.3I) . 

Second, if we use the classical Newton formula 

mMG 
F = = — , 



to calculate the galactic rotational velocity v r , then we have 
(7-5) f = 



Inserting (|7.4[) into (|7.5|l implies 

(7.6) u r = ==. 

V 1 + k r — kir 3 

As 1 > fco > fci, (|7.6j) can approximatively written as 

u r = V f 1 - ifc r + ifcgr 2 

which is consistent with the theoretic rotational curve as illustrated by Figure 
17.1( b). It implies that the distribution formula (|7.4|) can be used as a test for the 
revised gravitational held equations. 

Third, we now determine the constants fco and ki in (|7.2p . According to astro- 
nomic data, the average mass M ri and radius r% of galaxies is about 

M n = 10 n M ra ~2x 10 41 kg, 

(7 7) 

r x = 10 4 - 10 5 pc ~ 10 18 km, 

where M Q is the mass of the Sun. 
Taking V = 300km/s, then we have 

V 2 

(7.8) — = 8 x 10 23 kg/km. 
Based on physical considerations, 

(7.9) fc > hri (n as in ([777)1 ) 
By ([777 ]) - fT9 ]) . we deduce from that 

(7.10) fc - ^r^- - =4x 10- 18 km- 2 . 

Now we consider the constant k\. Due to the accelerating expansion of galaxies, 
the interaction force between two clusters of galaxies is repelling, i.e. for (|7.2[) , 

F > 0, r>f, 

where f is the average distance between two galactic clusters. It is estimated that 

f = 10 8 pc ~ 10 20 - 10 21 km. 

We take 

(7.11) f=-^xl0 20 km 

v2 
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as the distance at which F = 0. Namely, 



fc 1 



Hence we derive from (|7.10p and (|7.1ip that 

fci = k f- 2 = 10~ 57 km -3 . 

Thus, the constants fco and fci are estimated by 

(7.12) fc = 4 x 10~ 18 km~ 2 , fci = 10~ 57 km" 3 . 

In summary, for the formula (|7.2p with (|7.12[) , if the matter distribution M r is 
in the form 

(7 " 13) M -^TT^' 

then the Rubin law holds true. In particular, the mass M generated by the revised 
gravitation is 

V 2 V 2 n 

M = M T -M ri = —n-——^ , n as in HUD, 

G G 1 + k ri 

where Mt = V 2 r\jG is the theoretic value of total mass. Hence 

M fc n 4 

Mr ~ 1 + k n ~ 5' 

Namely, the revised gravitational mass M is four times of the visible matter M ri = 
Mt — M. Thus, it gives an alternative explanation for the dark matter. 

In addition, the formula (I7.2[) with (|7.12[) also shows that for a central held 
with mass M, an object at r > f (r as in (|7.1ip ) will be exerted a repelling force, 
resulting the acceleration of expanding galaxies at r > f . 

Thus the new gravitational formula (|7.2p provides a unified explanation of dark 
matter and dark energy. 

7.3. Effects of non-homogeneity. In this section, we prove that if the matter is 
homogeneously distributed in the universe, then the scalar potential <p is a constant, 
and consequently the scalar potential energy density is identically zero: $ = 0. 
It is known that the metric for an isotropic and homogeneous universe is given 

by 

dr2 +r 2 (d9 2 + sin 2 Odif 2 ) 



(7.14) ds z = -c z dt z + a 



1 — kr 2 



which is called the Friedmann-Lemaitre-Robertson- Walker (FLRW) metric, where 
the scale factor a = a(t) represents the cosmological radius, and fc takes one of the 
three numbers: —1,0,1. 

For the FLRW metric (|7.14l) . the energy-momentum tensor {T^} is given by 

(7-15) Tij = di&g(pc 2 ,g 11 p,g 2 2P,g33P), 

where p is the mass density, p is pressure, and 

Qj 2 2 2 2*2 /1 

.911 = 1 _ kr2 , 922 = a r , g 33 = a r sin 0. 
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By (|7.14|) and (|7.15|1 . the Einstein field equations are reduced to two equations: 

4ttG / 3p x 
- [ P+ V 2 



(7.16) a" = -^ [p+%)a, 



(7.17) aa" + 2{a'f + 2kc 2 = 4ttG (p - J) a 2 , 

and the conservation law div = gives 

,~ „<, 3 dR f p\ 

Then only two of the above three equations (|7.16p - (|7.18p are independent, and are 
called the Friedman equations. 

On the other hand, for the metric (|7.14j) with ()7. 15|) , the new gravitational field 
equations (|1.4j) with scalar potential are reduced to 

(7.19) a" = —\o+-£\a+-<p»a--a'<p' 



„ 20) '1 +2 (sL)\^ =1wG{i ,.^.^.^, 

a V. a / a \ c z / 2 2 a 

and the conservation equation 

a' ,, „ ( , a! a! p 

7.21 w'" + 3— w" = 87rG p' + 3— p + 3— 4 

a \ a a c z 

Again two of the above three equations (|7.19p - (|7.2ip arc independent. It follows 
from (f7~19)) and (|7T20| that 

(7.22) (a') 2 = ^-a 2 p - \^"a 2 - kc 2 . 

In the following we shall prove that 

(7.23) Vtp = 0, 

Namely, (p ^constant. In fact, let <p" and p have the form: 

(7.24) ,= 4 „" = 4. 

a J a a 

Inserting (j7T24|) in ((7~19l) . (f7~2T|) and (j7T22l) . we arrive at 

,,„ r , „ 4ttG IV 47TG,, 1 , , 

(7.25) a" = j + 7. ~ n°V ' 



(7.26) (or _^ + ^ = _^ > 

3 a 3 a 

(7.27) V = 8ttG0' + 247rGa 2 a'p/c 2 . 
Multiplying both sides of (|7.25l) by a' we obtain 

<« 8 » il 



. ,, 2 8nG9 lip 

(a ') + o- 

3 a 3 a 



4ttG(9' IV 4ttG , 1 



6 a c 2 " 2 



paa' - ^(a')V- 



It follows then from ([7T26| - ([7T28]) that 

1 
2 

which implies that (|7.23p holds true 



^(«')V = o, 
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The conclusion (|7.23p indicates that if the universe is in the homogeneous state, 

4 

then the scalar potential energy density gfg^ is identically zero: $ = 0. This fact 
again demonstrates that tp characterizes the non-uniform distribution of matter in 
the universe. 



8. Conclusions 

We have discovered new gravitational field equations (|1.4j) with scalar poten- 
tial under the postulate that the energy momentum tensor Ty- needs not to be 
divergence-free due to the presence of dark energy and dark matter: 

1 87rG7 

With the new field equations, we have obtained the following physical conclu- 
sions: 

First, gravitation is now described by the Riemannian metric <?jj, the scalar 
potential tp and their interactions, unified by the new gravitational field equations 



Second, associated with the scalar potential tp is the scalar potential energy 

4 

density jf-^^, which represents a new type of energy/force caused by the non- 
uniform distribution of matter in the universe. This scalar potential energy density 
varies as the galaxies move and matter of the universe redistributes. Like gravity, 
it affects every part of the universe as a field. 

4 

This scalar potential energy density -^q & consists of both positive and negative 
energies. The negative part of this potential energy density produces attraction, 
and the positive part produces repelling force. Also, this scalar energy density is 
conserved with mean zero: 

/ $dM = 0. 
Jm 

Third, using the new field equations, for a spherically symmetric central field 
with mass M and radius ro, the force exerted on an object of mass m at distance 
r is given by (see (|6.33|) ^l : 



F = mMG 



1 1 / 5\ , Rr 



o \ r I o 



R = $ for r > ro- 



where <5 = 2MG/c 2 . 

Fourth, the sum e = e\ + £2 of this new potential energy density 

1 8ttG 

and the coupling energy between the energy-momentum tensor Tjj and the scalar 
potential field tp 

c 4 /2 2MG\ dtp 
£2 ~ ~8^G \r + ^ r J ~dr' 
gives rise to a new unified theory for dark matter and dark energy: The negative 
part of e represents the dark matter, which produces attraction, and the positive 
part represents the dark energy, which drives the acceleration of expanding galaxies. 
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Fifth, the scalar curvature R of space-time obeys: 

R= 8 ^T + $. 

Consequently, when there is no normal matter present (with T = 0), the curvature 
R of space-time is balanced by R = $. Therefore, there is no real vacuum in the 
universe. 
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